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ABSTRACT 

Aims. The Coma cluster is the richest and most compact of the nearby clusters, yet there is growing evidence that its formation is still 
on-going. A sensitive probe of this evolution is the dynamics of intracluster stars, which are unbound from galaxies while the cluster 
forms, according to cosmological simulations. 

Methods. With a new multi-slit imaging spectroscopy technique pioneered at the 8.2 m Subaru telescope and FOCAS, we have 
detected and measured the line-of-sight velocities of 37 intracluster planetary nebulae associated with the diffuse stellar population of 
stars in the Coma cluster core, at 100 Mpc distance. 

Results. We detect clear velocity substructures within a 6 arcmin diameter field. A substructure is present at ~ 5000 km s"' , probably 
from in-fall of a galaxy group, while the main intracluster stellar component is centered around ~ 6500 km s" ' , ~ 700 km s" ' offset 
from the nearby cD galaxy NGC 4874. The kinematics and morphology of the intracluster stars show that the cluster core is in a 
highly dynamically evolving state. In combination with galaxy redshift and X-ray data this argues strongly that the cluster is currently 
in the midst of a subcluster merger, where the NGC 4874 subcluster core may still be self-bound, while the NGC 4889 subcluster core 
has probably dissolved. The NGC 4889 subcluster is likely to have fallen into Coma from the eastern A2199 filament, in a direction 
nearly in the plane of the sky, meeting the NGC 4874 subcluster arriving from the west. The two inner subcluster cores are presently 
beyond their first and second close passage, during which the elongated distribution of diffuse light has been created. We predict the 
kinematic signature expected in this scenario, and argue that the extended western X-ray arc recently discovered traces the arc shock 
generated by the collision between the two subcluster gas halos. Any preexisting cooling core region would have been heated by the 
subcluster collision. 

Key words. (ISM:) planetary nebulae: general; galaxies: cluster: general; galaxies: cluster: individual (Coma cluster); galaxies: 
evolution 

1. Introduction to fiy - 28mag arcsec^^, have shown that the ICL is made 

up of a wealth of diffuse features ranging from extended low 
Diffiise intracluster hght (ICL)__has iiow been pbs^erved surface brightness envelopes around giant ellipticals to long, 
m nearby ( |Feldmeier et al. | |2004; Mihos et al. 2005) thin streamers, as well as smaller scale features associated with 
and in intermediate redshift clusters (Zibetti etiiL 2005; ^^^^ yirgo galaxies. A similarly deep wide-field image of 
i Krick, Bernstein, & Pimbblet| |2006|). Individual intracluster j^e ICL in the Coma cluster is not yet available, but exist- 
stars have been detected in the Virgo and Coma clusters 



stars have been detected in the Virgo and Coma clusters j^g observations indicate similar complexity o n a variety of 

([Ferguson et al.| ll998t |Arnaboldi et al. | |2003t [Gerhard et al.| scales, including tidal features like streamers (iGregg & West! 

120051). Recent studies show that the intracluster hght contams arcs (&entham & Mobash eJ [Toll). extended h alos 

ofthe order of 10% of the mass in stars overall ( |Aguerriet al. | around D galaxies (,Baum et aI.„1986tlAdami et and 

| 2005MZibetti et al. || 2005|) , but m cores of dense and rich clusters extended background hght (iMelnick. White & Hoessel 19771: 



lIOOllZiMti^tiDSQOi, but in cores of dense and rich clusters extended background hg ht (IMelnick. White & Hoessel 19771: 

hke Coma, the local ICL fraction can be as high as 40-50% | Thuan & Kormendvl 19771) 
dThuan & Kormendvlll977l;lBernstein et al.lll995b . 

The large scale structure of the ICL in nearby clusters is Morphological studies of the ICL have some limitations 

rather complex. The recent surfac e brightness i neasur ements though: they cannot establish whether the ICL is made up by dif- 

of the ICL in the Virgo cluster by iMihos et akl (|2005|) . down fuse halos superposed along the LOS, but still physically bound 

to the galaxies they surround, or rather by stars that are free- 

* Based on data collected with the FOCAS spectrograph at the Aying i" the cluster potential. This question must be answered 

Subaru Telescope, which is operated by the National Astronomical by spectroscopic observations of the hne-of-sight velocities of 

Observatory of Japan, during observing run S04A-024. intracluster stars. Presently, the only tracers which allow us to 
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measure the kinematics and dynamics of the ICL are Intracluster 2. The velocity distribution of intracluster Stars in a 
planetary nebulae (ICPNe), whose [OIII]5007A emission can be Coma COre field from MSIS Observations 

used both for their identification and radial velocity measure- 

2.1. The MSIS technique 



By measuring the projected phase space of ICPNe we can 
constrain the dynamical age of the ICL compone nt, how and 
when this light origin ated (Napolitano et al. 2003; M urante et alj 
I2004t IWillman et alj|2004t iMurante et al.ll2007 ). Owing to the 
small fluxes of distant ICPNe this was so far only possible for the 
nearby Virgo cluster dArnaboldi et al.ll2004l) . With a new multi- 
slit imaging spectroscopy (MSIS) technique, essentially a spec- 
troscopic blind search technique, we have now been able to mea- 
sure ICPN velocities at substantially fainter fluxes. In our first 
application of this technique we have measured the velocity dis- 
tribution of ICPNe in a field in the inner core of the Coma cluster, 
approximately at the peak of the X-ray emission dGerhard et alj 
[2005;.Arnaboldi et al...2007,) . 



The Coma cluster (A1656) is the richest and most com- 
pact of the nearby clusters, and has been the subject of 
extensive study. Originally thought of as the prototype of a 
rich, regular, and relaxed cluster, it is now known to contain 
significant small-scal e and large-scale substr u cture in the 
galax y distribution (.Fitchett & WebsteJ Il987t iMelUer et al] 



119881 iBiviano et all 



. 1996h. galaxy vel ocity distribution 

dCofless & Dunn"1996^ 'Adami et alj|2()05bh. and X-ray emis- 
sivity (Briel, Henrv & Bohring^ 19921: IWhite. Briel & Henrvl 
119931: INeumann etalJI200U i2003h . Based on the recent work 
Coma is believed to be undergoing several accretion and merger 
events. Notably, the subcluster around NGC 4839 can be well 
distinguished in both the galaxy distribution and velocities and 
in the X-ray maps, and is falling into the main Coma cluster 
from the SW, at a velocity of ~ ITOOkms ' dMellier et alj 



119881: IColless&Duiinl 119961: [Neumann et all |200 lb . But also 
the cD galaxy NGC 4889 together with an associated subclus- 
ter are believed to be merging with the main Coma cluster, 
often associated with N GC 48 7 4, although the detai ls are less 
clear dC ofless & Dunn| Il996t lAdami et all l2005bl) . Finally, 
I Adami et aL C2005b) list a number of substructures, many of 
which may be infalling into the Coma cluster as well. 



The purpose of this paper is to ask what can be learnt by com- 
bining the measurement of the ICL kinematics in our field in the 
Coma core with the properties of the galaxy velocity distribution 
and X-ray morphology. Section |2] contains a brief summary of 
the MSIS technique and discusses the criteria for identifying the 
different emission sources in the data. We distinguish between 
galaxy PNe and ICPNe in the field, and discuss the velocity his- 
togram of ICPNe. The detailed procedures used for the data re- 
duction, the catalogue of PNe, and their spatial and magnitude- 
velocity distributions are given in Arna boldi et alj (1200 7). In 
Sect. [3] we discuss the velocity distribution of the ICPNe sample 
in relation with the galaxy redshift distribution and some notable 
features in the X-ray emission in the Coma cluster core, and de- 
velop a model for the on-going merger of the NGC 4889 sub- 
cluster with the second main subcluster in Coma around NGC 
4874. Finally, we give our conclusions in Section]?] In what fol- 
lows we assume a distance to the Coma cluster of 95 Mpc, thus 
1" = 460pc. 



Beyond about 20 - 30 Mpc distance, PNe are too faint to be de- 
tectable with narrow band surveys or slitless spectroscopy - their 
emission disappears in the sky noise in the narrow band filter. 
The brightest PNe in the Coma cluster at 100 Mpc distance have 
line fluxes of 2.2 x 10"'** erg s"' cm"^, corresponding to ~ 20 
photons per minute through an 8m telescope aperture, of which 
~ 2 will reach the detector for a typical ~ 10% overall system 
efficiency. To detect such distant and faint PNe requires a spec- 
troscopic blind search technique: spectroscopic, so that only the 
sky noise within a few A dilutes the emission from the PN, and 
blind, because the positions of these faint PNe cannot be previ- 
ously determined. 

The Multi-Slit Imaging Spectroscopy (MSIS) technique is 
such a spectroscopic blind search technique. Its first application 
to the detection of PNe in the Coma cluster was described by 
iGerhard et al.l d2005 ). The technique combines a mask of paral- 
lel multiple slits with a narrow-band filter, centered on the red- 
shifted [OIII]/15007A emission line. Spectra are obtained of all 
PNe that happen to lie behind the slits. The narrow band filter 
limits the length of the spectra on the CCD so that many slits 
can be simultaneously exposed. For each set of mask exposures 
only a fraction of the field is surveyed: to increase the sky cov- 
erage the mask can be stepped on the sky. The technique is sim- 
ilar to the approach us ed in the search for Lya emitting galaxies 
at very high redshifts dStern & Spinrad] 119991: ]TYan et al.ll2004t 
iMartin & Sawickill2004 . 

Our pilot survey in the Coma cluster was carried out with the 
FOCAS spectrograph at the Subaru telescope. The instrument 
was used with a mask of 70 parallel slits and a narrow band fil- 
ter of FWHM 60 A, centered on the [OIII] 5007A line at the 
redshift of the Coma cluster. The fraction of the field surveyed 
by a single mask is then ~ 12%. Each spectrum extends over 
43 pixels on the CCD, and the spectral resolution is 440km s 
As these values show, a compromise must be found in the MSIS 
setup between the number of slits and hence number of objects 
detected, and the velocity resolution and signal-to-noise (S/N) of 
the emission sources, which are dependent on the spectral reso- 
lution. For a detailed description of the ob serving technique an d 
the signal-to-noise calculations we refer to lGerhard et al.ld2005h . 

Three mask configurations were observed for the Coma 
field centered at Qf(72000) 12 : 59 : 41.784: 5(72000)27 : 
53 : 25.388, near the centroid of the X-ray emission in the 
Co ma cluster core , nearl y coincident with the field observed 
by Bernstein et al.l dl99 5') and ~ 5 arcmin away from the cD 
galaxy NGC 4874. The data reduction of the three masks ob- 
served with the MSIS technique was carried out in IRAF and 
is described in lArnaboldi et al.l d2007l) . The removal of the in- 
strumental signature was done following standard procedures for 
CCD reductions. The final dispersed image of the field resembles 
a brick wall made up o f adjacent 60 A wide two-dimensional 
spectra: see Fig. 2 o f iGerhard et alj d2005l) and Fig. 1 of 
lArnaboldi et al. I d2007l) . In this image point-like, monochromatic 
emission sources appear unresolved in both the spatial and wave- 
length directions, stars appear as spatially unresolved 60 A con- 
tinuum spectra, and galaxies appear as extended blobs. 
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2.2. Emission sources in tine piiot study MSIS fieid 

In the final co-added frames we look for emission line objects 
and classify them accordin g to the following criteria, described 
in mo re detail in .Gerhard et aD (12005.) and Arnaboldi et al.i 
(I2007h : 

- unresolved (both in wavelength and in space) emission line 
objects with no continuum, which are our intracluster plane- 
tary nebulae (ICPNe) candidates; 

- spatially unresolved continuum sources with unresolved or 
resolved line emissions. These are most likely background 
galaxies; 

- unresolved line emissions associated with the extended con- 
tinuum halos of Coma cluster galaxies in the field - such 
objects are compatible with being PNe associated with the 
stellar population emitting the continuum light. 

In total 60 such emission line objects were detected, of which 
according to these criteria 

- 35 FN candidates are ICPNe, with no detectable continuum 
flux down to 1.6 X 10"^"ergs"' cm^^A"', and equivalent 
widths ranging from EW >110Atol8A; 

- 20 are background objects; 

- 5 are PN candidates possibly associated with extended con- 
tinuum halos from nearby Coma galaxies. 

The identification of the objects unresolved in both position 
and velocities as PNe is further supported by the fact that the 
emission fluxes of the brightest objects are consistent with those 
of the brightest PNe in a population at distance 100 Mpc, and 
by the fact that both lines of the [OIII] doublet were seen in 
all four sources with sufficient flux > 4.0 x 10'^ erg s"' cm"^ in 
the [OIII] 5007 A line and at sufficiently large recession veloc- 
ity yobs > 7400 kms ' that also /14959A is redshifted into the 
wavelength range probed tGerhard et al.l 120051 ; I Arnaboldi et alj 
I2007h . Furthermore, as we shall see below, their distribution of 
recession velocities is centered around the Coma cluster and is 
inconsistent with a population of background objects uniformly 
distributed in velocity. It is possible that the sample contains a re- 
maining contamination from so far undetected, unresolved low- 
luminosity background emission line galaxies, but the contin- 
uum limit rules out compact HII regions such as or brighter than 
observed in Virgo jGerhard et al.ll20d2l) . The fi nal catalogue of 
PN can didates from the whole dataset is given in lArnaboldi et alj 
(I2007h . 

2.3. PNe in Coma Gaiaxies 

For the PN candidates superposed on the extended continuum 
halos of Coma galaxies, as well as for the ICPN candidates near 
such halos, it is necessary to check the measured LOS veloc- 
ity against that of the galaxy, in order to confirm or rule out 
the PN candidate-galaxy association. In the following, the IDs 
of the ICPN candidate s refer to their entries in the catalogue in 
lArnaboldi et"an(l2007l) . 

- The measured LOS velocities for IPN217 and IPN224, 6440 
and 6364 kms"' respectively, are consistent with the sys- 
temic velocity of the galaxy NGC 4876 onto which they are 
superposed, for which Vj, , = 6678 km s"'. 

- The measured LOS velocity of candidate IPN114 is 
6410kms"', not consistent with the systemic velocity of 
flie underlying galaxy CGCG160-235, v,v, = 8 114 kms"'. 



Similarly, for the two candidates IPN225 ( Vlos = 
5861 kms"') and IPN113 (vlos = 6719kms-'), the com- 
parison with the systemic velocities of the nearby Coma 
galaxies [GMP83-3376]v„., = 6815 kms"' and [GMP83 
3383]Vjv.s - 4640 kms"' shows very large difi'erences, of the 
order of 1000 kms"'. These three PN candidates are there- 
fore not bound to the galaxies, similarly as the Virgo ICPNe 
discovered serendipito usly along the line-of-sight to M86 
( lArnaboldi et al.lll996h . 
- However, close to IPN114 and CGCG160-235, one of the 
IPN candidates, IPN104, has Vlos = 8520 kms"' which, 
given the spectral resolution, is similar enough to the velocity 
of CGCG 160-235 that IPN 104 may be bound to CGCG160- 
235, even though it is not directly superposed on the galaxy's 
continuum image. Halos extended out to 70kpc in virial equi- 
librium have been detected in Vir go ellipti cals via surface 
brightness measurements (Mihos et alj2005h and ICPN Vlos 
measurements (lArnaboldi et al.ii2004h . 

In what follows, we consider the PN candidates IPNl 14, IPN225 
and IPNl 13 as intracluster PN, while IPN217, IPN224, IPN 104 
are considered as PN bound to Coma Galaxies. Thus the total 
number of ICPNe detected in our Coma field is 37. 

2.4. Veiocity distribution of ICPNe 

The LOS velocity distribution of all ICPN candidates is shown 
in Figure[T] this shows a clear association with the Coma cluster. 
The average velocity of the distribution is 6315 kms"' and the 
standard deviation is 867 kms"', but two substructures are also 
clearly visible: a main peak at Vlos ~ 6500 kms"' and a sec- 
ondary peak at Vlos ~ 5000 kms"'. The main peak in the LOS 
velocity distribution of PN candidates is very close to the sys- 
temic velocity of NGC 4889, V j,..? = 6495 km s"'. The secondary 
peak could be associated with a dissolved galax y from one of the 
subgroups discussed by lAdami et al.l (l2005bl) ; see Section 13.6! 
below. 

Fig. [U also shows the LOS velocity distribution of Coma 
cluster galaxies in a 6' x 6' field centered on the MSIS point- 
ing. Also this distribution has multiple peaks: the middle one is 
coincident with the v„,s = 6495 kms"' of NGC 4889, and two 
less prominent peaks are at Vlos < 5000 kms"' and Vlos ~ 
8 1 00 km s" ' . Thus the LOS velocity distributions of ICPNe and 
galaxies in the field more or less correspond; according to a 
Kolmogorov-Smirnov test, the probability of both being drawn 
from the same distribution is 16%. 

As a consistency check, we also show in Fig.[T]the distribu- 
tion of the observed emission lines of the background objects 
within the filter passband. For an easier comparison with the 
ICPN diagram, we express the Aobs in velocities, assuming that 
/}() = 5007 A. The distribution of point-like continuum sources 
with emission line is consistent with a flat distribution within the 
filter passband window used for the MSIS technique, and is con- 
sistent with a uniform background galaxy populatioifl There is 
no clear association with the Coma cluster here. The difference 
to the velocity distribution of Coma ICPNe is clear, and is one 
further argument in support of our source classification. 



' The apparent absence of a few objects on the blue edge of the filter 
bandpass is probably caused by a slight oflFset between the flat field 
image and the scientific image, which amplifies the noise on this side of 
the spectrum. 
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4000 5000 6000 7000 8000 9000 4000 5000 6000 7000 8000 9000 4000 5000 6000 7000 8000 9000 

Fig. 1. LOS velocities of ICPNe, Coma galaxies, and background emission galaxies in the MSIS field. Left: histogram of the LOS 
velocities of ICPNe in our Coma field. The thin red line shows spectra of a star as measured through the MSIS slits, giving the 
filter bandpass expressed in velocity. The velocities of the two supergiant galaxies NGC 4889 and NGC 4874 in the Coma core 
are indicated by vertical dotted lines; that of the nearer NGC 4874 is displaced from the main peak of the ICPN distribution. 
Center: LOS velocity histogram for the Coma galaxies in a 6' x 6' field centered on the MSIS field, whose velocities would fall in 
the velocity window allowed by the MSIS observing technique. The galaxy redshifts are from Adami et al. (2005b). Right: LOS 
velocity distribution of the line emitters associated with spatially unresolved continuum sources: these are compatible with a uniform 
population of background objects. Thin red line as in the left panel. 



3. Kinematics of diffuse light and the subcluster 
merger in the Coma cluster 

How does our measurement of the intracluster light kinematics 
in the Coma cluster core impact on the dynamical understanding 
of this densest and richest of the nearby clusters? 

Optical and X-ray data have shown that the Coma cluster 
contains major substructures around the infalling giant galaxy 
NGC 4839 as well as around the central pair of supergiant 
galaxies NGC 4874 and NGC 4889 (iFitchett & Webste3 fl987t 
ICoUess & Dumill 19961; iNeumann et all l2003h . In addition, a va- 
riety of smaller su bstructures have b een identified in position- 
radial velocity data dCoUess & Dunnll9 96; Adami et al. 2005 b). 
These findings, combined with the modern understanding of 
structure formation in the hierarchical Universe, have lead to 
the view that the cluster continues to accrete groups from the 
surrounding large-scale structure, and possibly has recently had 
a larger accretion of a subcluster centered around NGC 4889 
. The other supergiant galaxy in the core, NGC 4874, is often 
regarded as the true nucleus of the Coma cluster in these scenar- 
ios. This is based on the fact that NGC 4874 is surrounded by 
the largest concen tration in the projected distribution of galax- 
ies dBiviano et al.lll99 6). as well as to some extent on the fact 
that it has a relatively large diffuse light halo. The fact that its 
radial velocity is significantly offset relative to the velocities of 
nearby galaxies is then attributed to the interaction with NGC 
4889 , which could have per turbed NGC 4874 out of the cluster 
core (ICoUess & DunnllT996l) . 



3. 1 . The intracluster light 

How does the velocity distribution of ICPNe measured in our 
MSIS field some ~ 130kpc from NGC 4874 fit into this picture? 
And how can we use it to derive some more specific constraints 
on the merger history of the Coma cluster? The main peak of 
this velocity distribution is centered not at the systemic velocity 
of NGC 4874 at v = 7224 km s"', but around v = 6500 km s"', 
some 700 km s"' off and exactly at the systemic velocity of the 
other, more distant supergiant galaxy NGC 4889. Clearly, these 
ICPNe are not bound to NGC 4874 it self, the vel ocity dispersion 
of this cD galaxy is 284 km s"' (.Smith et al J200ft) . although they 




Fig. 2. The position of our MSIS field (red circle) on the diffiise 
light isodensity contours drawn by Thuan & Kormendv ( 197^ 
in the Coma cluster core. The MSIS field is about 5' south of 
NGC 4874. The second Coma cD galaxy NGC 4889 is 7' east 
(to the left) of NGC 4874. The bright object nor t h of N GC 4874 
is the star which prevented ' Thuan & Kormendvl ( 11977 *) from re- 
liably determining the northern parts of the isodensity curves in 
their photographic photometry. Note the strong elongation of the 
distribution of ICL in the Coma cluster core. The likely orbits 
of NGC 4889 and NGC 4874 up to their present positions are 
sketched by the yellow dotted and magenta dashed lines, respec- 
tively; see Section [ 



might be bound to the subcluster core around NGC 4874 if this 
is still present. 

Figure |2] shows the position of our MSIS field relative to 
the two s upergiant galaxies and the isodensity contours of ICL 
drawn bv lThuan & Kormendvl (Il977h . These early ICL measure- 
ments were photographic, but agreed well with previous pho - 
toelectric measurements by Me lnick. White & Hoessell ( Il977h . 
The contours are well-defined south of the two cluster dominant 



O. Gerhard et al.: Kinematic Substructures in the Coma Cluster Core 



5 



• I " ' I " " I "I" I " " I " "i 

NGC 4874 6'x6';^ ^ — 




Fig. 3. The two subcluster cores around NGC 4874 and NGC 
4889. Top-left: DSS image of the Coma cluster center. Bottom- 
left: positions of the galaxies in this region brighter than V=17; 
these are concentrated around NGC 4874 and NGC 4889. Right: 
velocity distributions for these bright galaxies that fall within a 
36 squai-e arcmin region around NGC 4874 (top) and NGC 4889 
(bottom). 



galaxies, while there are significant uncertainties on the north- 
ern side due to stray light from a bright star. In projection both 
galaxies appear to be embedded in a common envelope of dif- 
fuse light; as Thuan & Kormendy state, the light is not asso- 
ciated with any particular galaxy. The large elongation of this 
diffuse light distribution together with our kinematic measure- 
ment make a strong case that this ICL in the core of Coma is 
in a highly dynamically evolving state - we clearly do not see 
a relaxed ellipsoidal distribution of IC stars here. Moreover, the 
paucity of diffuse light stars in this field with NGC 4874 veloc- 
ities shows that the ICL in this area is poorly mixed. Mixing is 
expected to change the velocity distribution over a typical orbital 
time-scale, so the ICL in this area must have been released very 
recently. 

3.2. The subcluster cores 

Then the galaxies and subcluster cores from which the ICL is 
likely to derive will be in a stat e of str ong interaction as well. 
iKorchagin. Tsuchiva & Mivamal (1200 lb consider a head-on col- 
lision of the two cD galaxies and suggest that this could give rise 
to the elongated ICL distribution, as could a collision with mas- 
sive clust er substructures. Anothe r illustration is the early binary 
model of lValtonen & BvrdI (Il979i) . who estimated the masses of 
the subcluster cores around NGC 4874 and NGC 4889 and at- 
tempted to understand the dynamics of the Coma cluster in terms 
of the ef fects o f this binary on the surro unding galaxies. Later 
work bv Fitchett & Websteij (119871) and iBiviano etall d 19961) has 
confirmed the existence of subclumps in the galaxy distribu- 
tion centered on the two supergiant galaxies. The two dominant 
galaxies are projected onto the two mean peaks of the galaxy 
distribution. However, contrary to naive expectation, their veloc- 
ities do not mat ch the mean velocities of the peaks they are pro- 
jected onto. As IColless & DunnI (11996) state: "the positions of 
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Fig. 4. The Vlos distributions of Coma galaxies in 4 quadrants 
relative to the GMP83 Coma cluster center. Given in arcsec in- 
tervals along the RA and DEC directions, NE is [-600:0, 0:600], 
NW is [0:600, 0:600], SE is [-600:0,-600,0], SW is [0:600,- 
600:0]. In GMP83 coordinates NGC 4889 is at [-304,22], i.e., 
in the NE quadrant, and NGC 4874 is at [124,-41], i.e., in the 
SW quadrant. The MSIS field is in the SE, at [-15,-289]. 



the peaks and the dominant galaxies appear to be almost mirror- 
symmetric in velocity". 

However, this may be caused to considerable extent by the 
projection of foreground and background galaxies onto the sub- 
cluster cores. Such projection effects can be reduced by con- 
centrating only on the brightest galaxies, which are likely to be 
physi cally closer to the dominant galaxies NGC 4 874 and NGC 
4889 dValtonen & BvrdI 1 9791: iBiviano et aljl9"96l) . The distribu- 
tions and kinematics of the bright galaxies around the two super- 
giant galaxies are illustrated in F igure [3] which is based on the 
catalogue of Bivia no et al. I (11996). The top-left part of the fi gure 
shows a DSS image of the cluster center. In the bottom-left panel 
we show the positions of the galaxies brighter than V=17 in the 
region. Their concentration around NGC 4874 and NGC 4889 is 
apparent; a more quantitative analysis is in the references given 
above and will not be pursued here. The remaining two panels 
show the velocity distributions of these bright galaxies that fall 
in 6 arcmin squared regions around NGC 4874 and NGC 4889, 
respectively. The two histograms show an interesting difference: 
while the distribution of velocities around NGC 4874 is com- 
pact and may indicate a still self -bound cluster core, that around 
NGC 4889 is clearly not compact and has a twice larger veloc- 
ity width. The NGC 4889 subcluster core may thus currently be 
dissolving, consistent with the fact that the diffuse light in our 
MSIS field 8' west of NGC 4889 may also have been tidally dis- 
solved from the surroundings of NGC 4889, based on its similar 
radial velocity. 
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3.3. The NGC 4889 subcluster orbit 

Figure |4] shows the radial velo city histograms for all Coma 
galaxies from the catalogue of lAdami et alj (l2005bl) . in four 
quadrants of size (10 arcmin)^ NE, NW, SW, SE of t he cluster 
center defined by iGodwin. Metcalfe & Peac E (119831 hereafter 
GMP83). The strong variations from one histogram to the next 
show that the whole Coma cluster core is far from complete viri- 
alization. Note, for example, the strong peaks in the SE quadrant 
at vlos = 5000 km s"' and at 6500 km s"'; our MSIS field lies 
in the SE/SW quadrants (Fig. [3]). However, the real surprise of 
Fig.|4]is that it shows clear velocity peaks around the LOS veloc- 
ity of NGC 4889 in the SE and SW quadrants, with a further less 
significant peak in the NE quadrant, and clear peaks around the 
velocity of NGC 4874 in the NE and NW quadrants. A similar 
diagram showing the velocity histograms in (20 arcmin)^ sized 
quadrants shows a clear peak in the NE and a very strong peak 
in the SW around the NGC 4889 velocity, of which the SW peak 
is also robust against excluding the galaxies in the cluster cen- 
ter. Also, this diagram shows robust peaks around the NGC 4874 
systemic velocity in the NE and NW. 

This is contrary to expectation because in a general sub- 
cluster merger, these velocity distributions should have been 
changed and shifted by the gravitational pull during the inter- 
action. The only reasonable explanation for the histograms in 
Fig. |4] is that the average orbital accelerations on the galaxies 
of the two interacting subclusters on scales of ~ 500 kpc have 
predominantly been in the plane of the sky. Obviously this is a 
lucky coincidence; for a general orbit, the line-of-sight compo- 
nent of the acceleration would be visible in the distribution of 
galaxy recession velocities. For example, in a head-on collision 
along the line-of-sight, the observed galaxy velocities at the time 
of strong interaction would have been modified most in the inter- 
acting cluster cores and least in the subcluster outskirts. Then we 
should see a concentration of galaxies around the systemic ve- 
locities of NGC 4889 and NGC 4874 only in the cluster center, 
but not on larger scales. 

The idea that the two subclusters are meeting on an or- 
bit mostly in the plane of the sky is supported by two fur- 
ther arguments. First, the LOS velocity difference between the 
two supergiant galaxies and their surrounding cluster cores is 
only ^ 700 km s"', while the orbital collision velocity must ex- 
ceed ~ 2000 km s"'. To estimate this we have computed the 
virial velocity of the Coma cluster V200 - 1500 km s ~ ' fro m 
the virial masses and radii determined by iGeller et"an (119991) : 
iLokas & Mamonl (2003). The collision velocity when the cluster 
cores meet will be significantly larger than the median subhalo 
velocity found by Havashi & White (2006) from simulations, 
I.IV200, because this is determined for a distribution of subhalo 
orbital phases. If we take the relative velocity reached by 10% 
of their subhalos, this gives 1.55v2oo - 230 km s~'. Second, 
the C oma cluster is part of the "Great Wall" dGeller & Huchral 
11989 ). a large-scale structure of clusters and filaments that ex- 
tends approximately perpendicularly to the line-of-sight. In this 
structure Coma is connected through filaments to neighbouring 
large clusters, notably A2199, A1367, and A779, at approximate 
position angl es from North t hrough East of 80°, 250°, and 275° 
(see Fig. 3 o fiAdami et alj2005b) . E.g., the subcluster around the 
cD galaxy NGC 4839 appears to be infalling towards the main 
Coma cluster from the filament connecting Coma to A1367, at 
an estimated angle with respect to the line-of-sig ht of 74° and ve- 
locity in the plane of the sky of ~ 1750 km s~' dColless & DunnI 
[19961: iNeumann et aLllMl . 



The direction of infall of the NGC 4889 subcluster similarly 
is likely to have been from one of the filaments that connect 
the Coma cluster with the gr eat wall - but from which of the 
three filaments did it come? iBiviano et alj d 19961) construct a 
smoothed map of mean recession velocity Vlos for the fainter 
galaxies in their Coma sample. These galaxies are presumably 
further from the main site of gravitational interaction near the 
cluster cores, and thus ar e better tracers of the original subclus- 
ter motion. Figs. 12, 13 in lBiviano et aLld 19961) show an interest- 
ing gradient in the smoothed 2D velocity field along a direction 
35° inclined with respect to the RA axis (i.e., along PA 35° resp. 
215°). The mean Vlos of the faint galaxies is Vlos - 6700 kms"' 
at 10' -20' SW of the center and increases smoothly to above 
Vlos - 7000 km s"' 10' -20' NE of the center. Together with 
Fig.|4]we interpret this gradient in the sense that the NGC 4889 
subcluster now dominates WSW of the center, while the NGC 
4874 subcluster is more prominent ENE of the center. Based on 
the elongated ICL morphology, strong interaction has already 
taken place between the cluster cores, so that they must have 
already passed through each other. Then the observed gradi- 
ent requires that the NGC 4889 subcluster has last moved from 
ENE (the direction of A2199) across the center towards WSW, 
whereas the NGC 4874 subcluster must have fallen in from the 
West (the direction of A1367 or perhaps A779). 



3.4. Comparison with X-ray data and cluster simulations 

This conclusion finds independent confirm ati on from the re- 
cent X -ray analysis of l Arnaud et alj d200ll) and [Neumann et akl 
d2003h . These authors find a significant residual X-ray emis- 
sion in their XMM-Newton data from a large arc-like region 
between the Coma cluster core and the NGC 4839 group, and 
a temperature enhancement on the side of this structure to- 
wards the cluster core. There is also a positive residual on 
the south-east side at a similar distance from the center, but 
with a cooler temperature: this may be associated with a group 



ter collisions dRoettiger. Loken & Burns 1997: Taki zaw; 



|cy clus- 
m99^ 



iRicker & SarazinI 1200 It iRitchie & Thomasi i2002i) have shown 
that the compression and shocks generated in the collision give 
rise to large changes in X-ray luminosity, temperature, and emis- 
sivity contour shapes. After the collision of the cores of two 
similar mass clusters a strong arc-shaped shock is driven into 
the outer parts of the cluster's X-ray emitting envelope. The 
morphology of these arcs is reminiscent o f the morphology of 
the large western X- ray residual in Coma (lArnaud et al.ll200l] : 
iNeumann e"tani2003l) . 

While simulations with idealized initial cluster models show 
strong, regular, and symmetric features, cluster mergers from 
cosmological initial condit ions result in mu ch more irregular 
morphologies (Rowley. Tho mas & Kavll2004b . The morphology 
of their cluster merger 13, about 300 Myr after the collision of 
the cluster cores, is not unlike that observed for the Coma cluster. 
At this time, the emissivity contours are elongated along the di- 
rection of approach and the main arc-shaped shock is in the for- 
ward direction. From Fig. 3 of Rowley. Thomas & Kav (2004]), 
we estimate that this merger increased the mass of the main clus- 
ter by ~ 50%. Thus a picture in which the NGC 4889 subcluster 
has recently fallen in from ENE and its centroid is now past the 
center of Coma towards WSW, appears consistent with the pres- 
ence of the large western arc in the X-ray emission and the in- 
creased X-ray temperature associated with it. For reference, the 
western arc shock would have reached its current radius of ~ 30 
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arcmin from the cluster center after 300 Myr with an average 
velocity of ~ 2700 km s"'. 

We me ntion that our in t erpret ation is different from the sug- 
gestion by iNeumann et a n (I2003f) that the western arc could be 
caused by the hot plasma of a subcluster that previously merged 
into Coma from the southwest, which was heated in the colli- 
sion and deflected northwards by the main cluster hot athmo- 
sphere. In the hydrodynamic simulations, such deflections oc- 
cur when the hot athmospheres collide before the two cluster 
cores have met. This is unlike the situation in the Coma clus- 
ter where, as seen from the southwest, the two subcluster cores 
around NGC 4874 and NGC 4889 are ahead of the western arc. 
On the other hand, the NGC 4839 subcluster to the southwest 
has not yet reached t he main cluster: the X-ray observations 
jNeumann et alj l2003l) show a depression in the emission be- 
tween the NGC 4839 group and main Coma cluster (see also 
INeumann etaLll2001h . 

Our interpretation is also different from the scenario of 
IWatanabe et alJ (Il999l) . Based on an earlier, somewhat different 
ASCA temperature map, they suggested that the NGC 4874 sub- 
cluster had merged into the cluster from the southeast. There is 
no filament in this direction connecting to a neighbouring clus- 
ter. Also the elongation direction of the ICL distribution and the 
velocity differences of the two cD galaxies with respect to their 
surrounding galaxy concentrations cannot be explained if both 
galaxies have just passed eachother along this direction. 

Cluster mergers can disrupt cooling c ores if enough gas fr om 
the secondary cluster can reach the core (■ Gomez et al]|2002l) . as 
hap pens particularly in head-on collisions of equal-mass clus- 
ters dRitchie & Th omas 2002). With the near-head-on geometry 
of the collision indicated by the ICL and galaxy velocity mea- 
surements, it is likely that a preexisting cooling core in the clus- 
ter would have been heated by the collision, consistent with the 
absence of a cooling core in the Coma cluster now. 

3.5. Second passage of the subcluster cores and predictions 
for diffuse iight kinematics 

In cosmological simulations of structure formation, substruc- 
tures usually collide and merge along highly radial orbits. In a 
nearly head-on collision of two extended structures, the outer 
envelopes of both structures merge in a slow oscillatory fashion 
along an essentially unchanging orbital direction. By contrast, 
the interaction of the cores is faster, stronger, and may involve 
a significant impact parameter. In this way the cluster cores can 
have deflected each other strongly while the moti on of the outer 
clusters still looks apparently unperturbed (e.g., iMurante et al.l 
l2007h . 

This general picture is the likely explanation for the fact that 
NGC 4889 is apparently behind (towards the ENE) of the cen- 
troid of its subcluster along the direction of merging as deduced 
above. The elongated distribution of diffuse light requires that 
the galaxy has already undergone strong interaction in the core. 
Also, the IC star velocities show that significant diffuse light 
was lost west of its current position. All this is best explained 
if the current direction of motion of NGC 4889 is towards the 
east. I.e., the galaxy has gone around the other (NGC 4874) sub- 
cluster core in a highly eccentric orbit, moving presumably from 
north-east to south-west and then in a sharp turn towards south, 
east, and again north-east around NGC 4874. Such an orbit for 
NGC 4889 and a similar, symmetric orbit for NGC 4874 are 
sketched in Fig.|2] The diffuse light that we now see trailing the 
current position of NGC 4889 was probably lost during the or- 
bit around the NGC 4874 subcluster core between the first and 




Fig. 5. Spatial distribution of all cluster members with avail- 
able velocities. The X and Y axes denote RA and Dec, given in 
arcsec relative to the GMP83 Coma cluster center, with North 
up and East to the left. Green circles are Coma galaxies with 
6250km s"' < Vlos < 6850km s"', red triangles are galaxies 
with 7050 km s"' < Vlos < 8000 km s"', blue squares are galax- 
ies with 4800 kms"' < Vlos < 5350 kms"'. The positions of 
NGC 4874, NGC 4889 and the MSIS field are indicated on the 
plot. 

second close passages, as indicated on the figure. The material 
unbound from both galaxies shortly after their first close passage 
is likely to be the origin of the diffuse light lobes in the south- 
west and north-east. Note that little material is unbound from the 
colliding galaxies before the first pericentre passage; hence the 
absence of ICL with NGC 4874 velocities in our MSIS field. 

This predicts that in future observations of diffuse light kine- 
matics we should see NGC 4889 centered velocities west and 
south of NGC 4874 and towards the current position of NGC 
4889, and more NGC 4874 centered velocities towards the north- 
east of NGC 4874. The latter assumes that also NGC 4874 and 
its cluster core has already been stripped of significant light, 
which is plausible in the model above but less certain consid- 
ering the more intact nature of the NGC 4874 cluster core as 
discussed in Section IT2l 

One may ask whether a trace of the predicted ICL kinemat- 
ics around the two cD galaxies can be seen in the velocities of 
nearby galaxies dissolved from the cluster cores. Figure|5]shows 
the spatial distribution of galaxies with different velocities in the 
central parts of the cluster. In the distribution of galaxies with 
velocities +300 km s"' relative to the systemic velocity of NGC 
4889, there is a hint of a filament connecting this cD galaxy 
with the region south of NGC 4874, which is covered by the 
MSIS field. However, the noise expected in distributions with 
such small numbers of galaxies precludes any more definitive 
conclusion. 

3.6. The 5000 km s ' peak 

In addition to the dominant peak around NGC 4889's systemic 
velocity, the line-of-sight velocity distribution of ICPNe in our 
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MSIS field contains a secondary (blue) peak at velocities around 
5000 km s ' . There are also three galaxies with velocities around 
4500 km s ' in the field (Fig.lT] left and center panels, respec- 
tively). What could be the origin of these fast-moving objects in 
the Coma cluster core? 

The list of galaxy associations in lAdami et al. I (I2005bh con- 
tains two groups with low velocities and within less than 1 Mpc 
of the MSIS field: G7 (an extended association 0.3-1.2 Mpc east 
of the center, with 8 members, mean velocity 56 14 km s~'), and 
G9 (a compact group of 3 members 0.5 Mpc south-west of the 
center, with mean velocity 5710 kms~'). Relative to the mean 
cluster velocity of ~ 6850km s"' (Colless & Dunn 1996), these 
galaxies have velocities of ~ 1200km s"' towards the observer, 
while the ICPNe and galaxies in the blue peak of the MSIS 
field velocity distribution have even larger velocities towards us, 
~ 1 600-2700 km s-i. 

It is conceivable that the blue-shifted ICPNe and galaxies in 
the MSIS field and some or all of the galaxies in G9 and espe- 
cially those in G7 are part of a single structure falling through 
the Coma cluster from behind. The 8 blue ICPNe would be part 
of the remains of a galaxy shredded by the strong tidal field 
in the dense cluster core, corresponding to a total stellar mass 
of a few times 10^ Mq. This conversion assumes an a param- 
eter for the f irst magnitude in t he PN luminosity function of 
~ 10"^ (e.g.. Aguerri et al]|2005 ). and a mass-to-light ratio of 
~ 5. Because the MSIS field is located in the dense cluster core, 
the fast-moving objects in this field are probably located in a 
much deeper part of the Coma potential well than the galaxies 
in the G7 and G9 associations. This would explain their signifi- 
cantly faster line-of-sight velocities relative to the Coma cluster 
mean. To explain a present spread of this structure over / Mpc in 
the plane of the sky would also require differences in the velocity 
component perpendicular to the line-of-sight, of ~ 1 000/ km s ' 
over the past Gyr. Thus if all the galaxies of G7 and G9 were part 
of this structure, it would have had to dissolve early-on during 
its infall into the cluster. 



4. Discussion and Conclusions 

The measured offset of ~ 700 km s~' of the ICP Ne associated 
with the diffuse light in the iBernstein et al.l (1 19951) field, relative 
to the systemic velocity of the nearer giant elliptical (cD) galaxy 
NGC 4874, indicates that at least part of the diffuse light halo 
around NGC 4874 is not physically bound to this galaxy. 

From the substructure in the galaxy velocities, the disturbed 
X-ray emission, and the radio halo in the cluster, it is gener- 
ally belie ved that the Coma cluster i s in a state of on-going 
merging JWhite. Briel & Henrvl 11 9931: IColless & Duiinl Il996t 
iNeumann et al.ll2003t lAdami et a l. 2005b, and others). Here we 
have shown that also the LOS velocity distribution of the ICPNe 
in the center of the cluster suggests an ongoing merger of two 
substructures in the Coma core, which is not yet virialised. The 
two substructures are associated with NGC 4889 and NGC 4874. 

Starting with the distribution and kinematics of the diffuse 
light in the cluster core, we have developed a model for the ge- 
ometry of this subcluster merger that is based also on the distri- 
bution of Coma galaxy velocities and the presence of structure in 
the X-ray emissivity and temperature distribution. In this model, 
the NGC 4889 subcluster fell into the region of the Coma cluster 
from ENE, the direction of the filament connecting Coma with 
A2199. Whereas the infall direction of the NGC 4874 subclus- 
ter was from the west, probably the direction of the filament to- 
wards A1367. The two subcluster cores continue to disrupt each 



other through their mutual gravitational interaction, with the dis- 
ruption of the NGC 4889 subcluster core apparently in a more 
advanced stage. The inclusion of morphology and kinematics of 
the ICL in this model was vital because these data constrain the 
dynamical evolution in the cluster core itself. The ICL kinemat- 
ics predicted by the model in further fields can be tested with 
future observations, as described in Section [33] 

Currently, the two subcluster cores are past their second, 
close passage, during which much of the ICL seen through 
ICPNe must have been dissolved. This follows from the absence 
of any sign of mixing in the ICPN velocity distributio n, and 
strongly supports the conclusion of Muran te et al.l (|2007) that a 
large part of the ICL in clusters originates directly or indirectly 
from mergers of the most massive galaxies. 

Several Gyr in the future, the two supergiant galaxies and 
their subcluster cores will have merged. With the further arrival 
of the NGC 4839 subcluster from the SW, the cluster may not 
come to rest for significantly longer. The part of the galaxy dis- 
tribution centered around the NGC 4889 velocity in the 4 quad- 
rant diagram comprises a significant fraction of all galaxies in 
the Coma cluster core, perhaps 30%. This suggests that a truly 
major subcluster merger is currently taking place. Together with 
the other evidence for infalling substructures this suggests that 
Coma is forming now! 
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